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ABSTRACT 

Context. The far-infrared spectra of circumstellar envelopes around various oxygen-rich stars were observed using the 
ISO Long Wavelength Spectrometer (LWS). These have been shown to be spectrally rich, particularly in water lines, 
indicating a high H2O abundance. 

Aims. We have examined high signal-to-noise ISO LWS observations of the luminous supergiant star, VY CMa, with 
the aim of identifying all of the spectral lines. By paying particular attention to water lines, we aim to separate the 
lines due to other species, in particular, to prepare for forthcoming observations that will cover the same spectral range 
using Herschel PACS and at higher spectral resolution using Herschel HIFI and SOFIA. 

Methods. We have developed a fitting method to account for blended water lines using a simple weighting scheme to 
distribute the flux. We have used this fit to separate lines due to other species which cannot be assigned to water. We 
have applied this approach to several other stars which we compare with VY CMa 

Results. We present line fluxes for the unblended H2O and CO lines, and present detections of several possible v-z = 1 
vibrationally excited water lines. We also identify blended lines of OH, one unblended and several blended lines of NH3, 
and one possible detection of HsO" 1 ". 

Conclusions. The spectrum of VY CMa shows a detection of emission from virtually every water line up to 2000 
K above the ground state, as well as many additional higher energy and some vibrationally excited lines. A simple 
rotation diagram analysis shows large scatter (probably due to some optically thick lines). The fit gives a rotational 
temperature of 670l 2 gg K, and lower limit on the water column density of (7.0 ± 1.2) x 10 19 cm -2 . We estimate a CO 
column density ~100 times lower, showing that water is the dominant oxygen carrier. The other stars that we examined 
have similar rotation temperatures, but their H2O column densities are an order of magnitude lower (as are the mass 
loss rates). 

Key words, stars: AGB and post-AGB, stars: individual: VY Canis Majoris - supergiants - circumstellar matter 



1. Introduction The Infrared Space Observatory (ISO; iKessler et all 

■ Il996l ) satellite has revolutionised the study of circumstel- 

' Circumstellar envelopes form around red giant and super- lar envelopes by allowing access to all the spectral fea- 

giant stars during the later stages of their evolution. The tures in the mid- and far-infrared range. ISO spectra ob- 

mass loss from the stars results in an expanding envelope of served using the Sho rt Wavelength Spectrometer (SWS; 

material consisting of dust and molecules whose absorption |de Graaw et al.lll996l) an d Long Wavelength Spectrometer 

and/ or emission can be observed from infrared through ra- (LWS: IClegg et alj Il996h have previously been analysed, 

dio wavelengths. In oxygen-rich stars with high mass loss showing that the majority of the detected features in the 

rates, the cooling of the envelope is believed to be dom- FIR range are due to the rotati onal transition s of water (W 

inated by rotational l ines from water (H 2 0) vapour (e.g. Hya S WS: lNeufeld et all (I1996D; W Hya LWS:lBarlow et alj 



Uusttanont et al.lll994D . which has a rich far infrared (FIR) (plio^ Jj/c^^ et all ([19990; VY CMa 

spectrum. The FIR range also contains the rotational tran- SWS: Ne ufeld et alj (|1999f )). 
sitions of other simple molecules and molecular ions. 

In this paper, we present a dataset containing the ISO 
LWS spectra of the red supergiant star VY CMa, and sev- 

* Based on observations with ISO, an ESA project with in- eral other evolved stars ' The basic Parameters of these 

struments funded by ESA Member States (especially the PI spectra, such as identifications of lines over the wavelength 

countries: France, Germany, the Netherlands and the United range 45-196 /jm and their fluxes are discussed here. We 

Kingdom) with the participation of ISAS and NASA. concentrate on the results from VY CMa, as it shows the 

** Current address: Space Science Department, Rutherford spectrum with the highest signal to noise ratio and the 

Appleton Laboratory, UK line identifications are typical of several other stars in the 
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Table 1. Characteristics of the stars used in this paper. 



Star name 


Distance 
(pc) 


Mass Loss Rate 
(M yr" 1 ) 


Ref. 


VY CMa 


1100 


0.5-1 xlO" 4 


(1), (2) 


IK Tau 


300 


4.9xl0 -6 


(3) 


TX Cam 


390 


2.6xl0" 6 


(3) 


RX Boo 


160 


2.2xl0" 7 


(4) 


IRC+10011 (WX Psc) 


830 


6xl0~ 6 


(5) 


R Cas 


107 


3.4xl0~ 7 


(6) 



(1) Reid fc Me nten (in prepara t ion); IChoil (I2006T): (2) 
IDecin et all (120061); f3)l01ivier et all (120011) 1(4) iTevssier et aU 
i|2006l ); (5) IDecin et al.l (|2007l ); (6) iTruong-Bach et al.l (|1999l ) 

sample. The SWS line s pectrum for th i s star has previ- 
ously been presented by iNeufeld et alj (|1999D . The con- 
tinuum emi ssion of the SWS an d LWS spectra has been 
modeled by lHarwit et al.l (|200l[ ) and the LWS line emis- 
sion has been briefly described bv lBarlowl (|1999f) . VY CMa 
is thought to be in a late stage of stellar evolution, and 
appears extremely bright at infrared wavelengths. It has a 
mass loss rate estimated to be of the order of 5 x 10~ 5 to 
1 x 1 0~ 4 M(7) yr" 1 with a period of higher mass loss in the 
past (iDecin et al.ll2006|) . It is lo cated at the edge of a molec- 
ular cloud (jLada fc Reidlll978l ). and has a distance d erived 
from H2O and SiO masers of ~1.1 kpc from the Sun (|Choil 
2006, Reid & Menten, in preparation). In Sect. [3l we esti- 
mate the column densities of water and CO in the envelope 
of VY CMa to determine which is the dominant oxygen 
carrier. In Sect. [H we compare the results from VY CMa 
with several other stars that were also observed with the 
ISO LWS. The characteristics of these stars are shown in 
Table [1] 



2. Observations and data reduction 

A complete set of uniformly processed LWS grating spec- 
tra (using ISO observingmode L01) have been uploaded 
to the ISO D ata ArchivqB as hi ghly processed data prod- 
ucts (HPDP; iLlovd et alTl2003f ). This allows us to make 



a detailed comparison of consistently and accurately cali- 
brated spectral lines observed toward a sample of different 
stars. The objects, positions and ISO target dedicated time 
(TDT) numbers of the data we have used are shown in 
Tabled] 

To check the accuracy of the standard HPDP pipeline 
processing for our individual objects, we downloaded and 
manually reduced the LWS L01 spectrum of VY CMa, com- 
paring it with the HPDP data (to examine the effects of 
automatically removing glitches and dark current subtrac- 
tion). Both the test observations and the HPDP were ini- 
tially processed with the LWS off line pipeline (OLP) ver- 
sion 10. 

In order to reduce the data by hand, w e used the LWS 
Interactive Analysis (LIA; iLim et aTll2002l ) software to op- 
timise the dark current subtraction for each LWS detector. 
The absolute responsivity correction was not adjusted in- 
teractively because the OLP data already contains a bet- 
ter method of calculating the corrections to that available 
in the LIA software. Glitches were carefully removed from 
each scan by hand using the ISO Spectral Analysis Package 



see www.iso.esac.esa.int/ida/ 



(ISAP; [St urm et al.lll998| ). In addition, for ward and back- 
ward scans were compared and the regions where one direc- 
tion was affected by the response time of the detectors were 
removed (this can occur when the grating was scanned in 
the direction of increasing detector response from the edge 
of each band). A correction for fringing was applied using 
the LWS defringing algorithm within ISAP. The forward 
and reverse scans for each detector were then averaged. 

A careful comparison of the interactively reduced and 
HPDP spectra showed that there were only minor differ- 
ences in the final strength of the lines. This is mainly due to 
the defringing step which was only performed in the HPDP 
processing for detectors that the pipeline judged signifi- 
cantly affected (in the case of VY CMa, only for detector 
SW5, whereas interactively we judged defringing to make 
most difference to detector LW4). Further differences that 
may cause problems are that in the HPDP data reduction, 
both scan directions were simply averaged together (with- 
out accounting for the transient response of the detectors) 
and deglitching was carried out automatically. However, 
neither of these appear to have led to the introduction of 
spurious features in the HPDP spectrum. 

The largest difference between the two data reduction 
techniques is an improvement in the spectral shape of data 
from detector SW1 in the HPDP data. This is due to an 
additional correction in the HPDP pipeline that removes 
a double-peaked structure (with peaks at ~45 fim and 
~48 fiin) that is believed to be due to spuri ous features 
in th e SW1 relative spectral response function (jLlovd et alJ 
2003) . Using this correction allows us to have much greater 
confidence in the spectral shape observed using SW1 be- 
tween 43 and 50 fim. 

In conclusion, we determined that the HPDP pipeline 
did not introduce spurious artifacts that could be mistaken 
for lines, and that it had superior calibration for detector 
SW1. Therefore, we have used only the HPDP data for 
all of the stars in our sample. We applied multiplicative 
shifts to individual detector data to align their continua by 
examining the overlap region between pairs of detectors. 
However, this was not important for the spectrum of VY 
CMa because the absolute flux calibration applied by the 
HPDP already gave very good agreement between adja- 
cent detectors, with the only exception being SW1, which 
had an absolute flux ^20% higher than SW2. For the other 
stars, the shifts were more often required for the long wave- 
length detectors (up to ~30%). In order to examine the line 
fluxes, we subtracted the continuum level by fitting a 3rd 
order polynomial baseline to the data from each detector 
independently. 

Several well known spurious features that resemble lines 
remain in the spectra, and these are labeled in the plots 
presented here (e.g. features in absorption at 77 /jm using 
SW5 and 191 /zm using LW5, and features in emission at 
107 /im and 109 /jm using LW2; T. Grundy, private com- 
munication) . Some of these lines occur in the overlap region 
between two detectors and so can be ignored by selecting 
the optimum wavelength ranges to use. To provide a fur- 
ther check that the line features we have identified are real, 
we compared the spectra to observations of the asteroid 
Ceres, which should not contain any spectral line features 
(but does contain the spurious features). 

Where there were multiple LWS observations of the 
sample stars that showed good agreement, these were co- 
added to increase the signal to noise ratio. This applied 
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Table 2. Log of the L01 observations used in this paper. 
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TX Cam 


69501069 


3428 


10-10-97 


05 h 00 m 50.37 s , 56° 10' 52.5" 


RX Boo 


62300514 


1328 


31-07-97 


14 h 24 m 11.62 s , 25°42'l3.9" 




59701307 


1630 


05-07-97 


14 h 24 m 11.62 s , 25°42'l3.9" 


IRC+10011 


57701103 


2796 


15-06-97 


01 h 06 m 25.96 s , 12°35'53.2" 




57700513 


1910 


15-06-97 


01 h 06 m 25.96 s , 12°35'53.2" 


R Cas 


56801440 


2204 


06-06-97 


23 h 58 m 24.38 s , 51°23'l8.l" 




40 50 60 70 80 90 100 

Wavelength (pm) 



Fig. 1. LWS spectra (45-197 /im) of the stars used in this paper in units of X A f\. From top to bottom: VY CMa, IK Tau, 
TX Cam, RX Boo, IRC+10011 (xO.l) and R Cas (x0.05). The spectra of IRC+10011 and R Cas have been multiplied 
by factors of 0.1 and 0.05, respectively for clarity in the plot. The arrows indicate the true flux level for these stars before 
the multiplication factor was applied. Each spectrum was smoothed to reduce the noise before plotting. 



to IRC+10011 and RX Boo (see Table©. Although sev- 
eral other observations were also available for R Cas (TDTs 
26300712, 56801715 and 56801552), we have used only the 
one with longest integration time (TDT 56801440). 



The LWS beam size was ~80 , and the spectral reso- 
lution was 0.29 /im for detectors SW1- SW5 and 0.6 /im 
for detectors LW1-LW5 (|Grv et al.ll2003h . The RMS noise 
in the spectra varies for the different stars, but is roughly 
(10-20) xlO -20 W cm~ 2 pLmT 1 in the short wavelength de- 
tectors (measured between 87 and 89.5 /jm), and (0.5- 
2.2) x 10~ 20 W cm~ 2 /im -1 in the long wavelength detectors 
(measured between 140 and 142 /im). The wavelength bin 
width is 0.04 /im for the short wavelength detectors and 
0.13 /im for the long wavelength detectors. The final spec- 
tra before continuum subtraction are shown in Fig. [TJ plot- 



ted as A 4 f\. In this figure, the spectra have been smoothed 
with a window of 5 wavelength bins. 



3. Results: VY CMa 

In this paper we concentrate on identifying and character- 
ising the spectral lines in the LWS data between 45 and 
196 /im. The following sections present the line identifica- 
tions for VY CMa in this range. 



3.1. Water 

The SWS spectrum of VY CMa has alread y been shown 
to be rich in the rotational lines of water (Neufeld ct al. 
Il999l ). They identified at least 41 features due to water 
in the spectral range 29.5-45 /im, with transitions up to 
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Fig. 2. Continuum subtracted LWS spectrum of VY CMa with all water line transitions with upper levels up to 2000 K 
above the ground state marked below the spectrum (para: lower, ortho: upper of the two groups) and other lines marked 
individually above the spectrum. 



2939 K above the ground state. Many of the remaining 
weak features in the SWS spectrum can also be explained 
by including higher energy transitions of water. Here, we 
have analysed the LWS spectrum and carefully checked all 
of the prominent emission features. This shows that most of 
the detected lines can be attributed to rotational transitions 
in the vibrational ground state of water. 



The LWS spectra after continuum subtraction are 
shown in Fig. [5] The LWS spectral resolution is lower 
than achieved with the SWS, and many of the features are 
blended together. However, at longer wavelengths there are 
at least some gaps with few or no water lines in at all, giving 
the possibility to uniquely identify several other molecular 
species. The continuum level was estimated individually for 
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Fig. 3. An example of our water line fit for VY CMa for the region 59-68 /im. The final fit is shown by the solid black 
line, and individual water lines by the dashed line. Additional transitions above 2000 K in energy are marked. 



Table 3. Possible rotational lines from the vibrationally 
excited v 2 = \ state of water in the LWS spectrum of VY 
CMa. 



below 2000 K clearly show that higher energy transitions 
are present. These wavelengths are individually marked in 
Fig. (particularly 68-70 /im and 152-156 /im). 



Transition Wavelength Energy of Flux^ 

(/im) Upper State (K) (KT 19 W "cm" 2 ) 3.2. Fitting the water lines 



3o3 _ 2l2 

633-624 

4 2 3 "4l4 



182.36 
140.06 
115.71 



2491.8 
3284.2 
2745.2 



1.3 ±0.3 

1.2 ±0.2 

2.3 ±0.8 



a Errors do not include systematic uncertainty in the baseline 
subtraction (~30%) 



each detector by fitting a 3rd order polynomial through the 
lowest points in the spectrum as described in Sect. [2 For 
wavelengths greater than ~80 /im, there are sufficient gaps 
between emission features, and this method should give a 
reasonable estimate of the continuum. However, severe line 
blending occurs < 60 /im, and in this region the continuum 
level may have been overestimated because there is insuffi- 
ciently clear space between the lines. This means that the 
line fluxes may be systematically underestimated. However, 
the fluxes for lines at wavelengths < 60 /im are unlikely 
to be underestimated by more than a factor of a few, be- 
cause sufficient contamination of the continuum would re- 
quire very strong emission from lines lying >2000 K above 
the ground state and corresponding emission >80 /im is not 
seen. 

In Fig. [2] all the central wavelengths of all rotational 
ortho- and para- water transitions in its vibrational ground 
state are marked for energies up to 2000 K. Virtually ev- 
ery one of these transitions is matched by a feature in the 
data. This upper state energy is not a fundamental limit on 
the features detected in the spectrum, but adopted to ease 
of viewing the labels displayed at the LWS spectral resolu- 
tion. The inclusion of lines with higher upper state energy 
levels suggests that many of the low lying lines are prob- 
ably blended with other water features. Several spectral 
regions that are free of lines from transitions with energies 



In this section, we describe fitting the water lines in the vi- 
brational ground state. Due to the large number of blended 
Tines in the spectrum, any determination of the flux of in- 
dividual lines is likely to be rather uncertain. In addition, 
there will be systematic errors introduced in the line fluxes 
due to the uncertainty in the underlying continuum level 
at the shortest wavelengths. In order to estimate the line 
fluxes, we have used a method that simultaneously fits the 
water lines in the vibrational ground state up to a certain 
energy level (upper state energy less than 2000 K above 
ground) . The best fit is found by increasing the flux density 
at each line centre until the model reaches the data. The 
line wavelengths are fixed at their catalogue valu es (taken 
from the JPL line catalogue; iPickett et al.| [l998L and line 
widths fixed at the spectral resolution of the ISO detectors 
(0.29 /im for LWS detectors SW1-SW5 and 0.6 /im for de- 
tectors LWl-LW5; fGrv et al.|[2003h . Lines occuring within 
the wavelength range of each of the 10 LWS detectors are 
fitted together. 

The flux densities at the centres of all lines measured 
by the detector are increased by a step that is weighted 
(for each line) relative to the Einstein coefficient for spon- 
taneous ernissj^n^j^^^with values taken from the JPL cat- 
alogue; [PicketFer^nillll), and exp(—E u /kT weig ht), where 
E u /k is the energy of the upper state above ground in K, 
and T we i g ht is a single temperature applied to all lines. This 
ensures that lines with higher weight are given preference in 
blends. However, it does not necessarily mean that the final 
flux will be distributed following the weights - the actual 
shape of the data may allow a line with a lower weight to 
have a higher flux than one with a high weight. The weight- 
ing scheme does not treat ortho and para lines separately. 
The basic step size in the model is set to 0.01 times flux 
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density at the wavelength of the strongest line, and the 
individual step for all other lines are multiplied by their 
weights, 



I 



A-r 



odd 



(1) 



where I\ is the flux density at the centre of the line, w\ is 
the weight for that line, defined by, 



.4 



W\ 



id 



exp(-E u /kT weight ) 



MAX[A ul ] MAX[exp(-E u /kT weight )] 



and s is the basic step size, defined by, 
MAX[I A _ da4o ] 



100 



(2) 



(3) 



where the operation MAX finds the maximum value over 
the lines occuring within the wavelength range of the par- 
ticular LWS detector being fitted. 

A first estimate for the value of T we ight to use in this 
weighting function was determined from a rotation diagram 
(see Sect. I3.3[) using a fit weighted relative to A u i only, 



Aui 



MAX[A uZ 



(4) 



Figure [3] shows a fairly representative example of the 
final fit, for the region 59-68 /xm. This region clearly shows 
several higher energy lines that were not included in the fit. 

The advantage of this technique to estimate the line 
fluxes is that it is much faster than carrying out a multiple 
Gaussian \ 2 fit- Moreover, the flux in complicated blends 
is not overestimated (as it could be in a free Gaussian fit 
if the line shape is affected by other species). The disad- 
vantage is that it takes no account of the overall excitation 
of the molecule, except within blends where the lines are 
explicitly weighted according to A u i and T we ight ■ The final 
errors on each line flux are highly correlated in the various 
blended groups and therefore difficult to rigorously esti- 
mate. However, this fitting method provides a good starting 
point to investigate the spectrum, and represents the best 
that can be achieved using data observed with the available 
spectral resolution. 

In order to check the reliability of the automatic fit- 
ting routine, we compared the fit using several differ- 
ent weighting schemes: no weight, weighted relative to 
A u i only (equation [4} , and weighted relative to A u i and 
exp(—E u /kT we i g ht) (equation [2]). For the few completely 
unblended lines, we also carried out a manual fit using ISAP 
and in general this agrees within the errors with the auto- 
matic procedure. The unblended line fluxes for these dif- 
ferent fit methods are shown for VY CMa and IK Tau in 
Table |4] We estimate that the fitting errors for unblended 
lines are 15-20%. 

The final line fluxes determined from the fit for all lines 
up to 2000 K above the ground state (including blends) 
are shown in Tables IA.11 IA.2I and IA.3I for all of the stars 
investigated. In order to check the blended fits, we have 
compared our resu l ts for R Cas with those obtained by 
iTruong-Bach et all (|l999h . Our fitted values are m reason- 
able agreement (in a line by line comparison, many fluxes 
agree to better than 50%, and the remaining discrepen- 
cies can be explained by the fact that we have included 
more blends). We estimate that our blended line fluxes are 



probably of the correct order of magnitude (unless the as- 
sumptions used for the weighting scheme within the fit are 
wildly incorrect). 

We have also checked our fitted fluxes with the analysis 
of Macrc ker et al.l (|2008f ) who fitted water lines in a similar 



sample of stars. In general at long wavelengths where the 
lines are well separated, we find very good agreement in fit- 
ted fluxes. However, the difference in our analysis is that we 
have included the possibility of blends with higher energy 
lines. This is clearly necessary in VY CMa t o account for all 
of the detected features. In the analysis of iMaercker et all 
(2008]), they sometimes calculated model values lower than 
their quoted line fits to the data. This may be partly due to 
not identifying blends - some of our fitted fluxes which in- 
clude blends (where Maercker et al. did not include a blend) 
give values closer to their model predictions. On the other 
hand, in some cases, our simple treatment that includes 
many blended lines seems to underestimate the flux com- 
pared to their model, suggesting that we may have not dis- 
tributed the flux correctly. However, this comparison shows 
that in these stars, it is very important to consider that de- 
tected features may be due to blends of water transitions 
over a very wide range in energy, and it is not enough to 
merely consider the (strongest) lower level lines. 

3.3. Excitation of H2O 

In order to investigate the fitted line fluxes, we 
have used a simple rot ation diagram approach (e.g. 
iGoldsmith &; Langer| [l999). Assuming that all transitions 
are in local thermodynamic equilibrium (LTE), there is a 
linear relationship between the natural logarithm of the up- 
per state column density per statistical weight and the en- 
ergy above the ground state (i.e. following the Boltzmann 
distribution), 



In 



fju 



In 



Q 



T k 



(5) 



where N u is the column density in the upper state, g u is 
the statistical weight of the upper state, iVtot is the total 
column density, and Q is the partition function at tempera- 
ture T. A plot of \n(N u /g u ) versus E u /k can be fitted with 
a straight line with gradient 1/T and intercept ln(Ar tot /Q). 
We have calculated the column densities assuming optically 
thin emission lines using the relationship, 



N u = 



4ttA 
hcAuiQ 



l\—data 



ax 



(6) 



where fi is the solid angle subtended by the water-emitting 
region around the star (assumed to be a disk with diameter 
4") and J I\-data<iX is the measured flux. This choice of 
the source size is motivated by t he siz e of the CO emission 
region imaged by iMuller et alj (|2007ft . Our adoption of a 
uniform temperature across a region of this size is highly 
unrealistic and so our results should be taken as qualitative 
estimates (also see further caveats on the results below). 

Figure S] shows the rotation diagram for VY CMa. Even 
though we know that this model does not accurately de- 
scribe reality (VY CMa actually has a very complex envi- 
ronment - see Sect. [S]), and many of the lines are likely to 
be optically thick and/or subthermally excited, we think 
it provides an instructive starting point to view the fitted 
fluxes, particularly since the number of lines that we are 
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Table 4. Fitted fluxes for the unblended lines in the spectrum of VY CMa and IK Tau. 
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Transition 


Wavelength 


No Weight 


Weighted 


Manual 


IK Tau Weighted 




M 


Fit Flux 


Fit Flux 


Fit Flux_2 


Fit Flux 




(10~ 19 W cm" 2 ) 


(10" 19 W cm" 2 ) 


(10" 19 W cm" 2 ) 


(HT 19 W cm" 2 ) 


432 -321 


58.699 


15.5 


15.1 


27.8±5.4 


3.3 


6l6 _ 5o5 


82.031 


10.9 


11.0 


14.0±2 


2.5 


221-I10 


108.073 


7.2 


7.1 


9.2±1.1 


1.9 


432-423 


121.722 


3.1 


3.0 


3.8±1.0 


0.73 


4 2 3 "4l4 


132.408 


3.1 


3.0 


3.6±0.8 


0.49 


836-743 


133.549 


2.2 


2.1 


2.7±0.8 


0.54 


5l4-5o5 


134.935 


3.1 


3.0 


3.8±0.8 


0.58 


330-321 


136.496 


2.9 


2.8 


3.3±0.9 


0.73 


431-422 


146.923 


2.0 


2.0 


2.5±0.7 


0.29 


532-523 


160.510 


2.0 


1.9 


2.5±0.4 


0.55 


2l2~loi 


179.527 


5.3 


5.3 


6.3±0.7 


1.1 


221~2l2 


180.488 


2.6 


2.6 


3.4±0.4 


0.85 


633-542 


194.422 


1.5 


1.5 


1.8±0.3 


0.44 



a Errors do not include systematic uncertainty in the baseline subtraction (~30%) 
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Fig. 4. Rotation diagram resulting from fit to water lines 
in the LWS spectrum of VY CMa weighted by A u i and E u . 
Ortho water lines are shown as filled circles and para water 
lines as crosses. The best fit to lines with upper state ener- 
gies less than 2000 K (ignoring outlying points) is shown, 
giving a rotational temperature of 670ti3o K. 



Taking into account the above uncertainties, the fitted 
value for the total column density should be taken as a lower 
limit, and the temperature as a very rough estimate of the 
actual kinetic temperature (also bearing in mind that we 
expect a range of different thermal environments). The best 
fit straight line has a slope corresponding to a rotational 
temperature of T Iot =Q70tHo K ( the 

error was derived only 
taking into account the scatter in the data points). There 
does not appear to be a significant difference between the 
ortho- and para- lines. Due to the challenges described ear- 
lier regarding accurate flux estimation, we have discounted 
far outlying points by carrying out a 2a clip. These points, 
which can lie far above the best fit line, could be due to 
miss-fitting the flux of lines in blends. The rest of the scat- 
ter is due to real departures from the simple model as de- 
scribed above. More detailed modelling, beyond the scope 
of this paper is required to investigate these effects. The 
isotopic lines of H^O could be used to investigate the op- 
tical depths, but this will require future higher resolution 
observations (a useful limit on the H^O emission is difficult 
with the current LWS data). 

The best fit intercept point gives N to t/Q = (1.17±0.2) x 
10 17 cm -2 . The partition function a t 300 K is given in the 
JPL catalogue ()Pickett et al.lll998h as 178.12. The value 
scales as T 15 which leads to a partition function of 594.47 
at 670 K, and this gives a lower limit on the total water 
column density of (7.0 ± 1.2) x 10 19 cm" 2 . 



fitting is very large (141 lines used in the fit). The lower 
energy lines in Fig. |4] (< 1500 K) are unlikely to be com- 
pletely optically thin, and lMaercker et all (|2008f ) found that 
the water lines in their sample of stars were generally sub- 
thermally excited (although more thermalised in IK Tau 
than the other stars they looked at). 

Although most of the results cluster along a straight 
line, there is a large scatter. This is exactly the effect we 
should expect to see for a nonlinear molecule which has a 
mixture of optic ally thick and thin lines th at are not nec- 
essarily in LTE (|Goldsmith fe; Langerlll999l see their Figs. 
5 and 8). The lines least likely to be optically thick are 
those with upper state energies above 1500 K, although 
these lines are more likely to be blended, and so have large 
fitting errors. 



3.4. Vibrationally excited water 

In additi on to transitions from the vibrational ground state 
of water, Ncufcld et alj (|1999T ) also identified 4 features with 
water in its vibrationally excited V2—\ state. Rotational 
transitions of water in its ^2 = 1 state have also been ob- 
served at millimetre and submillimetre wavelengths, show- 
ing that the emission probably comes from a region clos e 
to the star (jMenten fc Melnickl [19891 : iMenten et aT1l2006l ). 
At the spectral resolution of the LWS grating, vibrationally 
excited water is much harder to separate. However, in the 
longer wavelength detectors there are three features that 
match with wavelengths of i/2=l rotational transitions. 
These are shown in Table |3l It is possible that other vi- 
brationally excited lines exist in the spectrum but all other 
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transitions with similar energies occur in blends. Also, the 
ground state transitions were identified first, which may 
mean that some vibrationally excited lines could have been 
misidentified as ground state ones (for example, the ^2=1 
2i2 — loi line occurs at 170.928 /j,m, close to the feature we 
have labelled as 155n — 16214 in Fig. [2]). Higher spectral res- 
olution observations would be required to unambiguously 
identify all of the vibrationally excited lines. 

iMenten et alj (|2006f ) detected two transitions of vibra- 
tionally excited water using the 12m Atacama Pathfinder 
Experiment telescope (APEX): 523 — 616 at 336.2 GHz and 
661 — 752 at 293.6 GHz. They used a model to calculate 
the optical depths and intensities of their measured lines, 
assuming that the transitions were thermalised. We have 
applied the same model to calculate the optical depth and 
predicted flux of our three detected FIR transitions. The 
line wavelengths and tr ansition parameters were taken from 
the JPL line catalogue (|Pickett et al.|[l998l ). If we make the 
same assumptions about the emitting region (T=1000 K, 
thermalised lines, uniform region within 0.05" from the 
star) and assume a line width of 20 km s _1 (the same as the 
APEX lines), we derive high optical depths of 500, 270 and 
600 for the three lines in Table [3] The predicted fluxes are 
a few 10~ 21 W cm~ 2 , which underestimate our measured 
LWS fluxes by factors of 36-86. 

The model pred ictions were closer for the sub- millimetre 
lines observed by IMenten et al.l |2006), underestimating 
their measured intensities by only a factor of 2.3 and 5, 
and a fact or of 3 for the 5so — 643 23 2.6 GHz transition ob- 
served by IMenten fc Melnickl (jl989). The underestimated 
intensity was explained by the fact that the 294 GHz line 
may be boosted by weak maser action. 

We ha ve also applied the m odel to the SWS lines ob- 
served by iNeufeld et alj (|1999D . predicting that the SWS 
lines should be stronger than the LWS lines, which is not ob - 
served. The SWS fluxes measured bv lNeufeld et alJ (|l999h 
are reproduced by the model within factors of 2-8. The op- 
tical depths for the SWS lines are predicted to be in the 
order of 10 4 . 

The model seems to be much worse at predicting the 
FIR fluxes than those in the sub-millimetre and mid- 
infrared. At the high optical depths predicted, the only pa- 
rameters in the model that affect the final fluxes are the 
temperature of the emitting region, the line width and as- 
sumed source size. The flux calculated by the model could 
be brought into better agreement by a combination of in- 
creasing the temperature and size of the emitting region, as 
the line fluxes are ~ linearly proportional to both. However, 
increasing the assumed emission region size has a large ef- 
fect on the predicted sub-mm main beam brightness tem- 
peratures. The excitation temperature does not have a large 
impact on the sub-mm lines, but a much higher tempera- 
ture would be needed than is ind icated by previous mea- 
surements ( see IMenten et al.ll2006l ). 

This either indicates that there is another problem with 
the simple model (for example, different excitation condi- 
tions for the FIR lines) or that the measured LWS fluxes 
are overestimated (possibly due to additional line blending 
which has not been taken into account here). Observations 
at higher spectral resolution are needed to confirm the 
fluxes, and more detailed modelling is required, taking into 
account both the sub-mm and infrared lines. 
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Fig. 5. Rotation diagram for the CO lines observed towards 
VY CMa. The dashed line is a fit to the 4 lowest energy 
lines using the manually fitted fluxes from Tabled 



3.5. CO lines 

In the LWS spectrum of VY CMa, the water line density 
decreases toward the longer wavelengths. CO is known to 
be an abundant molecule in O-rich circumstellar envelopes 
and many observations have been ma de via its millimetre 
and sub-millimetre transitions (e.g. iKemper et al.l [2003; 
iDecin et al.l |2006t IZiurvs et al1l2009D . The lowest energy 
CO transition in the ISO spectral range is J=14-13, with 
upper state energy of 580 K above the ground state. 

Several CO lines have previously b een reported in the 
LWS spectral ra nge toward R Cas (|Truong-Bach et all 
1999) and W Hya (jBarlow et al.lll996h . In addition gaseous 
CO absorption has been observ ed around 4.5 /im in the 
SWS spectra of O-rich stars (e.g. ISvlvester et al.lfl997T) . 

In the LWS spectrum of VY CMa, we observe some con- 
tribution to the spectrum from all CO lines from J— 14-13 
up to J=25-24. These lines are labeled in Fig.O The higher 
energy lines are difficult to identify due to blending with wa- 
ter lines and the best detection is of the lowest energy lines 
which are generally well separated. The fluxes of these lines 
at long wavelengths were determined in two ways and are 
shown in Table [5l As a first estimate, the CO fluxes were 
determined by fitting the excess emission after subtracting 
the water line fit described in Sect. 13.21 The flux was also 
determined by a manual Gaussian fit using ISAP, taking 
account of nearby blended lines with a multi-Gaussian fit. 
The automatic fit generally underestimates the flux com- 
pared to the manual fit, probably because water was given 
preference in the blends for the automatic fit (although the 
automatic fit is within the errors of the ISAP fit). The true 
flux is probably between the two values. 

We have carried out a rotation diagram analysis of the 
CO fluxes in the same way as described in Sect. 13.31 for 
H2O. This is shown in Fig. [5j with a fit to the 4 lowest 
energy lines (the fit ignores the fluxes for J = 19 — 18 and 
J = 21 — 20 which have very large error bars). This gives 
a rough estimate of the rotational temperature and total 
column density as we are only fitting a few lines. However, 
the results should be more reliable than for H2O, because 
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Table 5. Fitted fluxes for the CO lines observed above Table 6. Summary of other lines detected in the VY CMa 
120 /im for VY CMa. spectrum. 



Transition Wavelength Auto Manual ^^^^^^^^^^^^^^^^^ = ^^^^^^^^^^^^^^^^^ = 

(um) Fit Flux Fit Flux™ Transition Wavelength Comment 



(1(T 19 W cm" 2 ) (1(T 19 W cm" 2 ) „ (J£E± 

3-5^ i ^ , r\W ^TT " T—Q To" 1 /r> 1 CO 1 /I C 



J=14- 


-13 


185.999 


1.72 


1.96±0.4 


OH 2 n 1/2 J=3/2-l/2 


163.1/163.4 


blend with CO 


J=15- 


14 


173.631 


1.35 


1.88±1.3 


OH 2 n 3/2 J=5/2-3/2 


119.2/119.4 


well separated 


J=16- 


15 


162.812 


2.37 


2.69±1.4 


NH 3 K=2, J=3+-2" 


165.60 


well separated 


J=17- 


16 


153.267 


1.54 


2.00±0.6 


NH 3 K=2, J=3"-2 + 


169.97 


blend with H 2 


J=18- 


17 


144.784 


-Blend- 


-Blend- 


H 3 0+ K=l, J=2-2 


183.68 


possible identification 


J=19- 


18 


137.196 


1.06 


1.31±1.1 


2p 2p 
^ r 3/2~ r l/2 

U-line 


157.74 


blend with H 2 


J=20- 


19 


130.369 


-Blend- 


-Blend- 


52.5 




J=21- 


-20 


124.193 


1.60 


0.92±2.0 


U-line 


86.6 





a Errors do not include systematic uncertainty in the baseline 
subtraction (~30%) 



the high excitation CO lines are likely to be optically thin. 
We have used the manual ISAP fluxes from Table [5] and 
fitted a straight line weighted by their errors. This gives 
T Iot = 260 ± 140 K and N to t/Q = x 10 15 cm" 2 . The 

partition function a t 300 K is given in the JPL catalogue 
(|Pickett et"aT1ll998[ ) as 108.865, and this leads to a total 
CO column density of N tot = (9±? 4 ) x 10 17 cm" 2 . 

This column density is approximately 100 times lower 
than the H2O column density calculated in Sect. 13.31 show- 
ing that generally in the VY CMa envelope, water is the 
dominant oxygen carrier (although the H2O and CO come 
from different temperature regions). 

3.6. Other lines: OH, NH 3 and H 3 0+ 

Several other lines were detected in the spectrum of VY 
CMa. These are summarised in Table [75] and described in 
the following paragraphs. 

We detect strong emission from OH at 163 /im, which is 
blended with the CO J=16-15 line. This line is due to the 
lowest energy pure rotational transition in the 2 n 1 / 2 ladder 
of OH. Some contribution from the next transition in this 
ladder at 99 /im is probably also detected. In addition we 
detect emission in several other transitions of OH such as 
the lowest rotational line in the 2 n 3 / 2 ladder at 119 /im. As 
the cross ladder tra nsition at 34 11m is se en in absorption in 
the SWS spectrum (jNeufeld et al.|[l999] ). we might also ex- 
pect the other cross ladder transitions in the LWS range to 
be in absorption (e.g. the other transitions from the 2 n 3 / 2 
J=3/2 level at 53 /jm and 79 /im). However, these are both 
blended with water emission. 

The FIR lines of ammonia in its ^2 bending-inversion 
mode arc rotation-inversion transitions, with the lowest 
energy transition K=2, J=3-2 at 165.60 and 169.97 /im. 
There is a clear emission peak at 165.60 fim in the spectrum 
of VY CMa (see Fig. [2]). This does not overlap with any 
water lines and has a width comparable to the instrument 
resolution. We think that the assignment of this feature to 
NH 3 is secure because it does not appear to be blended. 
There are only a few other species that could provide alter- 
native explanations for the line: the most likely two choices 
are HNC J=20-19 and H 3 0+ K =0, J=5-4. The HNC line 
has a high energy upper level at 913 K and the correspond- 
ing transition in HCN at 161.35 /im is not clearly detected. 
If the line was due to H 3 + , there should have been other 



features present which are not observed. The line is also 
observed in the spectrum of IK Tau (see Sect. [4]). 

The NH 3 line at 169.97 /im occurs in a blend with water 
lines. The transition at the next highest energy in NH 3 is 
K='S, J=4-3 at 124.65 and 127.11 /mi. Both these lines 
occur in the wing of water emission (but are marked in 
Fig. E]). 

Ammonia has been observed in the envelopes of both 
carbon and oxygen rich stars - the first detect i on in an 
O-rich envelope was made by iMcLaren fc Bet j (|1980t ) in 
VY CMa, VX Sgr and IRC+10420. They observed several 
vibration-rotation transitions in the V2 band at 10.5 and 
10.7 /im. Subsequently these transi tions have been mea- 
sured in more detail for VY CMa by[M 

onmcr et~all (p2000h 

and references therein. These observations show that NH 3 
probably forms near the termination of the gas acceleration 
phase. 

The only detection of the radio inve rsion transitions of 
ammo nia in O-rich stars was made by [Menten fc AlcoleU 
(|1995l ) toward IK Tau and IRC+10420. They found that the 
strongest inversion transition for IK Tau was K=3, J=3- 
3 whereas for IRC+10420 was K=l, J =1-1. This implies 
that with a higher mass loss rate, photodissociation occurs 
further from the star and the NH 3 emission comes from a 
lower temperature region. 

The H 3 + ion has a similar pyramidal structure to 
the isoelectronic molecule NH 3 . In an analogous way to 
NH 3 , the oxygen atom can tunnel through the plane of 
the molecule leading to an inversion splitting of the ro- 
tational levels. However, for H 3 + , the inversion split- 
ting is very large (55.3 4 cm -1 for the z/2 bending-inversion 
mode; iLiu et all [l986). This means that the pure inver- 
sion transitions occur at FIR wavelengths rather than in 
the radio regime as for ammonia. Several rotation-inversion 
transitions at 1 mm h ave previously been detected in gi- 
ant m olecular clouds ( W ootten et al.l Il991l : [Phillips et alj 
1992 ) and the lines near 300 GHz mapped toward Sgr B2 
(|van der Tak et aTl |2006[ ) . Several of the FIR transitions 
have also been observed in absorption toward Sgr B2 with 
ISO (jGoicoechea fc Cernicharol 12001b iPolehampton et all 
[200l . 

The peak abundance of H 3 + in oxyge n-rich circum- 
stella r envelopes is predicted to be ~10~ 7 (Ma mon et alj 
119871 ), with the major formation route being the photodis- 
sociation of H2O and OH. Excitation is probably via the 
absorption of mid-IR radiation rather than collisions and 
thus absorption bands may also be observable at 10 and 
17 /im. The only (possible) detection of H 3 + toward an 
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evolved star so far is t he 2a feature observed toward VY 
CMa by Phil lips et al.l (|1992f ) at the correct velocity for the 
K=0 J=3-2 396 GHz transition. 

In the LWS grating spectrum of VY CMa, the spectral 
resolution is low and the inversion and rotation-inversion 
transitions are difficult to separate from the water lines. 
However, there is an unidentified feature at the correct 
wavelength for the K=l J=2-2 transition at 183.68 fxm. 
The other low energy inversion transition may contribute 
to the (water) emission features around 180-181 /im. The 
rotation-inversion transitions are harder to positively iden- 
tify as they occur at shorter wavelengths where the density 
of water lines is higher. 

In several O-rich star spectra, a feature at 157.8 /im 
has been assigned to the ato mic fine structure l i ne of 
ionised carbon at 1 57.74 /im (|Truong-Bach et ail Il999t 
ISvlvester etahl Il997t) . This line also appears in the VY 
CMa spectrum but no other atomic lines are clearly visible. 
ISvlvester et alJ (|1997| ) attribute this line to bad subtrac- 
tion of the galactic background in their sources. However, 
at least some contribution to it can be explained by the 
955-862 para water line at 157.88 /mi. The upper level for 
this transition occurs 2031 K above the ground state. In 
VY CMa there could be some contribution from both lines. 

There are several remaining unassigned features in the 
spectrum and these are labeled as 'U-lines' in Fig. [5] In 
addition, there are several species which have transitions 
in the range and are predicted to be reasonably abundant 
in O-rich circumstellar envelopes, but are not possible to 
detect due to line ble nding with water (e.g. HCN, H 2 S; 
IWillacv fc Mllia^[l997h . 

4. Comparison with other stars 

The spectrum of the oxygen rich Mira variable IK Tau 
shows remarkable similarity to that of VY CMa. However, 
the line fluxes shown in Fig. [6] are ~4 times weaker and 
the continuum level ~10 times lower. Since it was only ob- 
served for about twice the integration time of VY CMa, 
the signal-to-noise of the spectrum appears lower than that 
of VY CMa. Within the noise level and accuracy of the 
continuum subtraction, all of the lines observed in the VY 
CMa spectrum are present in the IK Tau spectrum, gen- 
erally with very similar relative intensities. Figure [6] shows 
the spectrum of IK Tau scaled by a factor of 4.5 and com- 
pared to VY CMa. It is clear that the same lines are de- 
tected in each object. There are, however, a few exceptions 
where there are marked differences: there is much weaker 
OH emission with respect to water in IK Tau. The OH line 
at 119 /im is not present and there may be less contribu- 
tion to the 163 /im CO/OH blend from OH. Also the line at 
99 /im is weaker in IK Tau. However, the other lines such as 
CO, NH3 and the possible v-i = 1 water lines are all present 
in both stars. In addition, the unidentified feature at 86 /im 
is present in both spectra. There is only one feature that is 
present in IK Tau but not in VY CMa at 73.8 /im. 

In Fig. [6l we also plot the other stars from Table Q] 
In order to compare the spectral lines, we have multiplied 
each continuum subtracted spectrum by factors of 4.5 for 
IK Tau and R Cas; 6 for TX Cam and 10 for RX Boo 
and IRC+10011. We have also smoothed the noisy spectra. 
Several key spectral lines are shown. 

The other stars of a similar type to IK Tau (TX Cam, 
RX Boo, IRC+10011) have a much lower signal to noise 



Table 7. Results of fitting a straight line to the rotation 
diagram for all of the stars in our sample. 



Star 


Trot 


iVtot(H 2 0) 
(10 18 cm" 2 ) 


Number of Fitted 




(K) 


Data Points 


VY CMa 


670I 130 


70 ± 12 


141 


IK Tau 


5901$ 


1.5 + 0.3 


140 


TX Cam 


5401* 


3.7 + 0.9 


125 


RX Boo 


5001* 
1400™ 


3.0 ± 0.6 


125 


IRC+10011 


5.6 + 1.3 


130 


R Cas 


6501?™ 


8.4 + 2.3 


132 



ratio and their lines are weaker than IK Tau. However, the 
majority of detections in VY CMa and IK Tau also seem 
to be present in these stars. 

We have fitted the water lines in the vibrational ground 
state for the other stars in the same way as for VY CMa 
and the results are shown in Tables IA.11 IA.2I and IA.31 We 
have carried out a rotation diagram analysis for each star as 
described in Sect. 13.31 and the results are shown in Table [71 
giving a rough estimate of the rotational temperature and 
a lower limit for the total water column density. 

The rotational temperatures for the other stars are very 
similar (within 25%) to the value found for VY CMa, which 
indicates that the H 2 0-emitting gas is excited under similar 
conditions. If the excitation is dominated by collisions, T ro t 
should reflect the kinetic temperature of the gas because of 
the high densities of AGB star envelopes. Alternatively, ra- 
diation may dominate the excitation, in which case T lot 
would reflect the temperature of the ambient radiation 
field. The only star where T rot deviates from the trend is 
IRC+10011, but in this case, the error on T rot is rather large 
(due to a larger scatter in the data points in the rotation 
diagram) . 

The derived H2O column densities for the other stars 
are lower than for VY CMa, with values ranging from 2% 
to 12% of the value for VY CMa. Since the other stars also 
have 10-100 times lower mass loss rates (see Table Q}, this 
result indicates that water is also the main oxygen carrier 
in the other star envelopes. However, a direct scaling of 
A r (H 2 0) with mass loss rate docs not appear: the ratio 
of H2O column density and mass loss rate shows a large 
variation over our sample of stars (factor of >10). 

5. Summary 

In this paper we present the ISO LWS spectrum of the lumi- 
nous supergiant star VY CMa and give a brief description 
of the detected lines. We then compared the spectrum with 
several other evolved star spectra, particularly the Mira 
variable IK Tau. 

The detected lines in the spectra of VY CMa and IK 
Tau are remarkably similar, both in the species detected 
and the relative intensity of the water lines. 

We draw the following conclusions: 

— We confirm that the spectra of our sample of evolved 
stars are dominated by water line emission in the FIR, 
and we report the detection of nearly all the H 2 lines 
up to ~ 2000 K above the ground state. 

— A simple rotation analysis of the water lines shows 
that there are a probably a mixture of optically thick 



E. T. Polehampton et al.: The ISO LWS line spectrum of VY CMa 



11 



Frequency (THz) 

5.0 4.9 4.8 4.7 4.6 4.5 4.4 4.3 4.2 4.1 4.0 3.9 3.8 3.7 3.6 3.5 3.4 



-VY CMo 

K Tau x 4.5 
TX Cam x 6 

RX Boo x 10 
RC1001 1 x 10 

R Cas x 4.5 




140 

Wavelength (/^m) 



Fig. 6. Comparison of the spectrum of VY CMa with the other stars (the spectra of the other stars have been smoothed 
and shifted for clarity - see text). 



and thin lines which may be subthermally excited (the 



fit gives a rotational temperature of 
a corresponding column density -/V to t 



670;^ K, and 



(7.0 ± 1.2) 



10 cm which should be treated a lower limit to the 
true value). 

We tentatively assign several of the detected features to 
excited water in its 2/2 = 1 vibrational state. 
We estimate the column density of CO to be (9^ 2 ) x 
10 17 cm -2 , showing that water is the dominant oxygen 
carrier in these envelopes. 

We present a detection of ammonia in VY CMa and IK 
Tau, and a tentative detection of the H30 + ion toward 
VY CMa. 



We note that the kinematics around VY CMa are 
very complex as sho wn by interferometric CO observations 
(|Muller et al.l l2007h which suggest the presence of spheri- 
cally symmetric outflow together with a high velocity bipo- 
lar flo w. The water maser observations by M enten et alj 
( 2008) further imply that some fraction of the water emis- 
sion may be emitted by special regions associated with 
shock fronts, excited in photochemical reactions associated 
with the outflow material. The presence of this complex 
excitation environment has also been indicated by obser- 
vatio ns of other molecules (jZiurvs et alj|2009l : ISmith et alj 
2009) . To bring together an accurate understanding of the 
local environment will require a sophisticated model be- 
yond the scope of the present study. Future observations 
with Herschel and ALMA should bring the necessary data 
to understand this complicated source. In particular, HIFI 
will add the ground state water lines, the isotopic H^O 
lines and the line profiles which may be used to study kine- 



matics and self- absorption effects, and PACS will reobserve 
all of the lines above 55 /im. ALMA will give the spatial 
distribution of H2O and other molecules in these AGB en- 
velopes. 
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Appendix A: Full list of fitted fluxes 



E. T. Polehampton et al.: The ISO LWS line spectrum of VY CMa 
Table A.l. All line fluxes calculated by our fitting procedure for the 6 stars. 



Transition Wavelength Energy of VY CMa IK Tau TX Cam RX Boo IRC+10011 R Cas 





H 


Upper 
State (K) 


(10" 16 W m- 2 ) 












SW1 


















660-735 


43.790 


1175.0 


146 


29.7 


7.2 


5.0 


18.1 


11.1 


541 -4 3 2 


43.893 


550.4 


49.2 


11.5 


1.3 


5.6 


3.6 


5.4 


743 ~7l6 


44.048 


1013.2 


208 


49.9 


7.9 


21.1 


21.7 


17.0 


542-431 


44.195 


552.3 


183 


48.7 


17.7 


6.2 


12.1 


25.2 


836-725 


44.702 


1125.7 


174 


46.1 


7.3 


9.7 


19.1 


7.6 


523~4l4 


45.112 


323.5 


323 


102 


18.8 


33.4 


19.0 


32.5 


331-2o2 


46.484 


100.8 


200 


52.3 


11.4 


13.8 


19.7 


27.4 


642-615 


46.539 


781.1 


55.3 


22.4 


2.2 


5.7 


5.0 


5.2 


633-606 


46.544 


642.7 


13.0 


4.4 


0.49 


1.1 


0.82 


1.1 


1029-918 


46.608 


1552.6 


13.2 


5.5 


3.4 


0.83 





5.6 


735-624 


46.746 


867.3 


102 


17.4 


10.1 


0.11 


2.4 


17.7 


4 4 o-4i3 


47.028 


396.4 


63.1 


14.1 


4.5 


6.3 


2.0 


4.8 


10i9-9 28 


47.029 


1554.4 


4.0 


1.0 


0.26 


0.47 


0.24 


0.30 


II111-IO010 


47.040 


1603.6 




















lloii-lOno 


47.045 


1603.6 


117 


30.0 


7.5 


14.9 


8.2 


9.4 


844-919 


47.267 


1324.0 








1.3 











541-5l4 


47.420 


574.7 


312 


62.4 


24.5 


27.0 


21.2 


43.8 


761-836 


47.601 


1447.6 


115 


16.1 


8.1 


16.4 


2.0 


17.2 


532-423 


47.973 


432.2 


304 


64.4 


23.7 


23.1 


11.2 


37.7 


440-331 


49.282 


410.4 


6.2 


0.12 


0.01 


4.1 


0.02 


0.02 


661-734 


49.334 


1212.0 


77.2 


25.8 


2.2 


6.8 


5.9 


7.2 


441-330 


49.337 


410.7 


50.3 


15.2 


1.6 


3.9 


2.3 


4.8 


634-523 


49.391 


642.4 


167 


50.7 


13.3 


9.4 


4.5 


15.7 


928-8l7 


50.634 


1270.3 


260 


77.3 


9.9 


10.1 


16.4 


13.8 


SW2 


















927-836 


51.071 


1447.6 


267 


62.1 


15.1 


17.5 


9.2 


12.8 


10no-9o9 


51.445 


1323.9 


217 


51.2 


0.62 


10.2 


10.9 


4.8 


10oio _ 9i9 


51.461 


1324.0 


0.20 


0.05 





1.0 


0.01 


1.7 


918-827 


51.685 


1274.2 


230 


41.0 


18.5 


8.4 


6.3 


3.1 


550-625 


52.864 


795.5 


44.1 


14.1 


13.7 


7.9 





17.2 


533-422 


53.138 


454.3 


175 


37.7 


23.3 


12.1 


3.4 


11.1 


743-8l8 


53.455 


1070.7 


24.9 


12.8 


29.3 


10.9 


0.70 


8.3 


532 -5o5 


54.507 


468.1 


158 


60.6 


21.6 


9.2 


16.0 


14.5 


827-7l6 


55.131 


1013.2 


149 


34.2 


0.62 


26.6 


12.6 


9.8 


IO29-IO110 


55.840 


1603.6 





0.40 





6.6 


2.0 





651-726 


55.858 


1021.0 


104 


42.0 


0.54 


12.9 


3.8 





IO19-IO0IO 


56.027 


1603.6 


66.8 


1.9 


1.5 


0.03 








431-322 


56.325 


296.8 


155 


27.4 





14.9 








9l9-8o8 


56.771 


1070.5 


126 


12.9 


0.78 


13.6 








9o9~8l8 


56.816 


1070.7 


13.4 


3.9 


0.90 


0.94 





0.78 


853-928 


56.972 


1554.4 


109 


20.4 


0.15 


11.1 


4.6 


12.4 


752-827 


57.394 


1274.2 


73.4 


9.2 


7.2 


6.5 





9.7 


422~3l3 


57.637 


204.7 


142 


29.4 


10.6 


21.6 


6.1 


14.6 


8l7~726 


57.709 


1021.0 


221 


45.8 


10.4 


16.4 


13.9 


19.0 


642-7l7 


58.377 


843.8 





0.20 





5.3 








432 -321 


58.699 


305.2 


151 


33.4 


7.6 


1.5 


11.7 


15.0 


SW3 


















726 _ 6l5 


59.987 


781.1 


76.4 


17.7 


3.8 


4.8 


8.0 


1.8 


826-735 


60.162 


1175.0 


106 


28.4 


4.8 


2.5 


7.5 


7.2 


762-835 


60.229 


1511.0 


8.4 


0.90 


0.08 





1.8 


0.18 


541-616 


61.316 


643.5 


10.5 














1.0 


4 3 i-4 4 


61.809 


319.5 


95.2 


24.8 


4.9 


4.5 


2.5 


11.7 


4 4 o-5i 5 


61.916 


469.9 


1.9 


0.49 





0.06 


0.05 


0.06 


936-845 


62.418 


1615.3 


10.7 


3.6 


0.46 


0.31 


0.28 


0.49 


928~9l9 


62.432 


1324.0 


103 


27.8 


2.0 


0.41 


2.9 


1.4 


9l8~9o9 


62.928 


1324.0 


68.6 


16.4 


1.8 





5.5 


6.2 


8l8 _ 7o7 


63.324 


843.5 


71.5 


13.1 


3.9 


2.9 





11.7 


8o 8 -7i7 


63.458 


843.8 


139 


35.7 


7.2 


11.6 


3.6 


13.7 


762-753 


63.880 


1524.6 


12.1 


2.7 


0.01 


0.01 


0.23 


0.19 


661-652 


63.914 


1278.5 


88.9 


25.3 


3.1 


6.2 


4.5 


5.0 


660 -651 


63.928 


1278.6 


5.5 


1.6 


0.20 


0.34 


0.30 


0.30 


761-752 


63.955 


1524.9 


110 


32.4 


4.3 


5.2 


6.0 


5.0 


625-5l4 


65.166 


574.7 


107 


24.2 


1.3 


8.3 


3.6 


10.1 


716-625 


66.093 


795.5 


60.3 


27.2 


2.7 


2.2 


0.84 


8.6 


330-221 


66.438 


194.1 


156 


37.8 


10.3 


11.8 


5.2 


20.7 


331-220 


67.089 


195.9 


87.7 


20.9 


11.0 


7.0 


7.7 


8.1 


330~3o3 


67.269 


196.8 


153 


41.7 


6.3 


10.3 


5.5 


18.5 
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Table A. 2. Line fluxes for all stars, continued. 



Transition 


Wavelength 
(pm) 


Energy of 

Upper 
State (K) 


VY CMa 
(10" 16 Wm- 2 ) 


IK Tau 


TX Cam 


RX Boo 


IRC+10011 


R Cas 


SW4 


















827-818 


70.703 


1070.7 


81.7 


7.5 


6.6 


3.6 


7.7 


2.2 


524-4l3 


71.067 


396.4 


145 


25.9 


3.2 


8.3 


3.0 


3.6 


7l7~6o6 


71.540 


642.7 


56.2 


5.8 


1.2 


3.4 








551-624 


71.788 


867.3 


1.10 


0.57 








0.09 


2.1 


7o7~6l6 


71.947 


643.5 


85.4 


30.4 


5.4 


2.2 


4.4 


10.6 


817-808 


72.032 


1070.5 


18.3 


10.3 





0.81 


6.2 





937-928 


73.613 


1554.4 


37.1 


15.8 





2.0 


6.6 


1.4 


725-634 


74.945 


933.7 


81.0 


25.6 


3.7 


5.3 


6.1 


12.7 


321~2l2 


75.381 


114.4 


94.8 


23.0 


10.9 


10.2 


2.3 


14.7 


854-845 


75.496 


1615.3 


68.2 


30.6 


6.6 


3.8 


6.5 


12.7 


551-542 


75.781 


877.8 


47.2 


13.0 


0.69 


6.7 


5.0 





753-744 


75.813 


1334.8 


18.7 


4.7 


0.16 


2.9 


1.7 





652-643 


75.830 


1088.8 


0.03 


0.01 














550-541 


75.910 


878.2 


93.5 


33.8 


2.5 


0.84 


5.3 


5.2 


651-642 


76.422 


1090.3 


50.9 


15.9 








3.9 


5.3 


752-743 


77.761 


1339.9 


37.4 


10.1 





4.2 


3.9 





SW5 


















423~3l2 


78.742 


249.4 


110 


23.2 


11.2 


4.0 


1.7 


12.6 


6l5~524 


78.928 


598.8 


125 


26.1 


4.8 


7.8 


4.7 


9.8 


946-937 


80.222 


1749.9 


32.3 


6.3 


2.9 


5.2 








853-844 


80.557 


1628.4 


28.7 


6.2 


1.1 











726~7l7 


81.216 


843.8 


34.0 


4.8 


2.0 





0.84 


0.32 


9 27 -9i8 


81.405 


1552.6 


43.2 


10.2 


1.7 





0.06 





835-744 


81.690 


1334.8 


27.6 


2.4 


0.09 


4.1 


2.0 





6l6 _ 5o5 


82.031 


468.1 


110 


24.8 


8.8 


9.4 


1.4 


16.0 


836-827 


82.977 


1274.2 


35.3 


9.0 


5.0 


4.6 


5.6 


6.8 


6o6-5i5 


83.284 


469.9 


60.0 


17.0 


4.9 


3.7 


0.50 


6.0 


7l6~7o7 


84.767 


843.5 


75.2 


16.9 


5.1 


6.7 


7.6 


3.8 


845-836 


85.769 


1447.6 


36.5 


3.4 





4.1 


5.1 





322-211 


89.988 


136.9 


85.4 


19.0 


8.1 


7.5 


2.9 


10.0 


744-735 


90.050 


1175.0 


50.1 


10.8 


2.0 


10.4 


5.3 


11.8 


LW1 


















643-634 


92.811 


933.7 


43.2 


5.2 





3.0 


3.5 





735-726 


93.383 


1021.0 


40.3 


4.6 


1.1 





4.2 


3.0 


652-725 


94.172 


1125.7 


32.8 


5.0 


1.7 





3.9 


1.8 


542-533 


94.210 


725.1 


3.3 


0.66 


1.2 





0.35 


0.20 


625-616 


94.644 


643.5 


43.4 


8.1 


1.3 


1.9 


2.3 


2.9 


4 4 i-4 3 2 


94.705 


550.4 


65.3 


13.7 


3.0 


3.7 


4.0 


7.3 


945-854 


95.176 


1805.9 


46.6 


9.3 


1.9 


3.7 


2.9 


1.8 


5l5~4o4 


95.627 


319.5 


73.7 


14.0 


3.4 


4.9 


4.5 


10.8 


4 4 o-4 3 i 


95.885 


552.3 


59.1 


13.6 


3.9 


2.3 


3.1 


2.8 


937-10010 


98.329 


1603.6 


45.8 


9.3 


2.1 


2.3 


3.7 


2.3 


541-532 


98.494 


732.1 


54.8 


7.7 


1.1 


0.56 


3.7 


3.2 


5o5~4l4 


99.493 


323.5 


85.2 


18.8 


6.4 


6.0 


3.9 


15.4 


826-817 


99.979 


1270.3 


92.4 


22.0 


2.8 


3.6 


3.9 


6.1 


5l4~423 


100.913 


432.2 


89.3 


23.6 


6.0 


5.5 


3.7 


11.2 


220-lll 


100.983 


53.4 


105 


25.4 


5.9 


7.2 


3.9 


14.3 


624-533 


101.209 


725.1 


84.4 


21.2 


5.4 


3.5 


2.8 


8.7 


642-633 


103.916 


951.8 


52.8 


15.7 


2.4 


4.3 


3.6 


5.9 


615-606 


103.940 


642.7 


0.56 


0.17 


0.03 


0.85 


0.03 


0.07 


634-625 


104.094 


795.5 


35.5 


8.1 


3.4 


0.85 


1.7 


3.7 
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Table A. 3. Line fluxes for all stars, continued. 



Transition 


Wavelength 
(pm) 


Energy of 

Upper 
State (K) 


VY CMa 
(10" 16 Wm- 2 ) 


IK Tau 


TX Cam 


RX Boo 


IRC+10011 


R Cas 


LW2 


















221-I10 


108.073 


61.0 


70.6 


18.2 


9.7 


6.3 


4.1 


15.3 


524~5l5 


111.628 


469.9 


30.0 


8.5 


3.6 


0.74 


2.2 


3.4 


743 _ 734 


112.511 


1212.0 


21.8 


5.1 


1.5 


1.9 


2.1 


4.0 


4 4 i-5i4 


112.803 


574.7 


12.7 


3.0 


0.89 


1.7 


1.6 


1.7 


4l4~3o3 


113.537 


196.8 


59.2 


15.7 


5.1 


3.6 


3.0 


13.3 


533-524 


113.948 


598.8 


59.1 


13.8 


5.5 


5.1 


4.1 


10.7 


927-lOlK) 


114.454 


1603.6 


11.7 


5.1 


0.76 


1.1 


2.8 


1.9 


836-909 


116.350 


1323.9 


10.2 


2.2 





0.04 


0.92 


0.19 


734-643 


116.779 


1088.8 


21.3 


6.0 


0.35 


2.0 


1.4 


4.8 


946-853 


117.684 


1807.0 


1.7 








0.19 


0.52 


0.36 


937-844 


118.405 


1628.4 


11.0 


3.1 


1.2 


1.3 


0.81 


10.7 


432-423 


121.722 


432.2 


30.0 


7.3 


1.7 


1.7 


2.8 


4.5 


844-835 


122.522 


1511.0 


10.2 


2.2 


0.27 





0.55 


1.8 


9 3 6-9 2 7 


123.460 


1729.3 


10.7 


2.0 


0.22 


1.5 


1.7 


1.8 


4o4~3l3 


125.354 


204.7 


51.8 


15.3 


4.3 


4.0 


4.0 


7.3 


LW3 


















331-322 


126.714 


296.8 


33.0 


8.9 


3.5 


0.76 


4.1 


4.4 


725"7l6 


127.884 


1013.2 


14.8 


4.0 


2.7 


0.78 


2.4 


1.5 


945-936 


129.339 


1845.8 


3.0 


0.38 


0.61 





1.4 





753-826 


130.319 


1414.2 


8.4 


2.3 


1.0 


0.85 


0.45 


0.41 


4 2 3"4l4 


132.408 


323.5 


30.1 


4.9 


2.5 


2.4 


1.8 


4.3 


836-743 


133.549 


1339.9 


21.0 


5.4 


2.0 


3.3 


1.6 


1.7 


5l4~5o5 


134.935 


468.1 


30.3 


5.8 


2.0 


2.7 


0.99 


3.5 


330-321 


136.496 


305.2 


27.6 


7.3 


2.8 


3.1 


1.6 


4.0 


735-808 


137.683 


1070.5 


6.6 


2.1 


0.92 


1.4 


1.2 


1.8 


3l3"2o2 


138.528 


100.8 


45.1 


9.3 


6.0 


5.4 


1.8 


8.1 


844-753 


138.641 


1524.6 


11.4 


2.4 


1.1 


1.0 


0.53 


2.1 


4l3 _ 322 


144.518 


296.8 


32.2 


7.4 


3.1 


2.4 


1.4 


3.7 


431-422 


146.923 


454.3 


19.7 


2.9 


2.3 


0.9 


1.3 


1.5 


LW4 


















835-826 


148.708 


1414.2 


5.6 


0.83 


0.65 


1.4 


0.73 


0.71 


542-615 


148.790 


781.1 


3.6 


1.4 


0.29 


0.22 


0.65 


0.69 


322~3l3 


156.194 


204.7 


5.7 


1.5 


0.75 


0.67 


0.21 


0.54 


523-432 


156.265 


550.4 


37.0 


10.0 


2.9 


1.7 


1.6 


3.4 


331~4o4 


158.312 


319.5 


13.4 


0.48 


1.6 


0.33 


0.19 


0.83 


845-752 


159.051 


1524.9 


10.4 


1.8 


0.81 


0.59 


0.73 


1.2 


634-707 


159.400 


843.5 


6.8 


1.2 


0.79 


0.20 


0.45 


0.49 


826~9l9 


159.485 


1324.0 


0.76 


0.25 


0.70 


0.01 


0.14 


0.09 


532-523 


160.510 


642.4 


18.7 


5.5 


2.6 


0.09 


1.2 


1.6 


734-725 


166.815 


1125.7 


5.2 


0.58 





0.06 


0.33 


0.17 


624-615 


167.035 


781.1 


10.9 


3.7 


0.80 


0.26 


1.0 


0.82 


LW5 


















735-642 


169.739 


1090.3 


15.5 


0.19 


0.46 





1.4 


4.1 


633-624 


170.139 


867.3 


35.4 


9.9 


2.2 


0.55 


2.1 


1.4 


533-606 


174.607 


642.7 


48.8 


14.4 


4.3 


4.9 


2.6 


6.7 


3o3 _ 2l2 


174.626 


114.4 


3.2 


1.4 


0.34 


0.25 


0.10 


0.71 


432~5o5 


174.920 


468.1 


5.7 


0.10 


1.6 


0.63 


2.3 


0.60 


2i2-loi 


179.527 


34.2 


52.7 


11.2 


5.0 


5.5 





8.8 


221 _ 2l2 


180.488 


114.4 


26.4 


8.5 


2.1 


0.26 


3.7 


4.2 


4l3-4 4 


187.111 


319.5 


15.1 


4.9 








0.84 


1.3 


854-927 


187.810 


1729.3 


3.7 


0.91 





1.8 


0.43 


3.6 


643-716 


190.437 


1013.2 


6.8 


0.64 


1.4 





2.3 


1.3 


633-542 


194.422 


877.8 


14.6 


4.4 


1.5 


4.8 


0.88 


7.9 



